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A study of the Camargue rice field responses from Landsat 
imagery is made, taking intc. account the ecological and agricultural 
conditions of the zone. 


The large dispersion of the rice-fields responses' in comparison 
with other types of ground objects, permits to develop a fast method of 
classification using the linear classifier. 

The results are compared with that obtained by a conventional 

method * 
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I - Description of rice parcels in the Camargue 

The arrangement of rice fields in the Camarque (Southern 
France) depends on the ecological and agricultural conditions of the 
zone. 

Areas of rice cultivation consist of adjacent parcels whose 
dimensions depend mainly on wind conditions and topography. 

The parcels are leveled so that the depth of the water can 
he controlled to within one centimenter. Prior to levellings a topogra- 
phical study of the terrain is made. The parcels are delineated in such 
a way to minimize the volume of earth to be removed and to reduce 
the harmful effect of soil rasping. This implies that the parcels be 
narrow and that their length follow concur lines. 

Furthermore, since the prevailing wind, the Mistral, which 
causes damage particularly during the first months of vegetation, blows 
North-South, the fields are narrow and orientated East-West in order to 
reduce the action of the wind. 

Given these restrictions, most of the fields have a small 
area of from 1 to 3 hectares and are rectangular. 

Figure 1 shows a view of the Camargue rice fields taken on 
June 20 1975 by Daedalus multispectral Scanner at 7000 m (near Infrared 
Channel : 0,900 - 1,100 pm) 
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Figure 1 

View of the Camargue rice fields 
Scale : 1/100.000 * 

The fields are separated by dirt roads which must permit 
transportation operations (sowing, chemical treatments, harvest...) 
which often take place when the fields have been flooded. These roads 
are from 4 to 10 meters wide and take up from 15 to 25 % of the total 
cultivated area. 
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After f looding s the fields look line a mosaic of rectangular 
surfaces of about 1 to 3 hectares, with the water level at about 10 cm. 

During sowing, the layer of water serves as a temperature 
control for the shoots. Thus the incessant motion of the water helps to 
compensate for variations in air temperature. 

From August to September, the movement of the water .(emptying 
and refilling) depends on evaporation, losses through drainage and the 
water need of the plants. 

Therefore, irrigation condition differs from a parcel to another. 

Furthermore, during the vegetative stage of rice, the coverage of 
the vegetation on parcels at a given time varies from parcel to parcel, 
depending on the variety planted. t 

II - Study of the response of reference rice-fields 

The reference rice fields located within a block of parcels 
by means of ground observations and aerial photography are selected from 
Landsat 2 recordings. 

Each pixel on Landsat image either coincides with one parcel 
or straddles two. 

“ The homogeneity of the rice class ‘is. studied by means of the 

histograms of the responses of the selected samples. 

A zone is defined as homogenous when the gradient at all points 
in the zone is zero or below a threshold determined ' by the measurement 
noise , 

The gradient can be calculated following 8 directions. This 
leads to consider zones of 3 x 3 points (or pixels) . 
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the zone of 3 x 3 points is considered to 
be homogenous. 

In addition to this concept of instrumental homogeneity there 
is also the concept of the homogeneity of the ground object class. The 
latter concept is related to problems of resolution. 

With LANDSAT resolution (80 m x 80 m) such agricultural terrain 
as vineyards or orchards (with a structure of alternating vines and soil 
or trees and soil) appears as stable composite and the samples taken 
within this type of terrain have homogenous responses. 

However, for rice-fields in the Camargue, which consist of 
small parcels with differing states, differences in irrigation conditions 
and in rice varieties, the responses of the pixels taken within a sample 
area give a greater-dispersion, especially in the MSS7 channel. 

Figure 2 gives 2 examples of histograms for 2 of the most 
important types of agricultural terrain in the Camargue, rice fields and 
vineyards 
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Histograms of rice-field and vineyard samples. 
(HSS7 responses - July 6 1975) 


6 - Method of classification 


The fastest method for mapping of rice fields from data 
obtained from May to July, when ricefields are flooded or in early 
vegetative stage consists of density slicing of the MSS7 image. 


When data are acquired during the vegetative stage of rice, 
i.e* from July to September, responses in MSS7 channel of rice— fields 
overlap that of other crops (com, cereals), as it is shown, in the 
figure 3 , 


In the figure 3, three classes of objects are studied by a 
bi spectral diagram of the responses of samples in channels MSS5 and MSS7 


rice fields, vineyards and cereals {com, wheat) from the July 6 1975 
scene (2165-09511). 

The average variances in responses 'and the correlation coefficients 
•for the 3 elasses are the following : 



s 

<3“ 7 

r 

Rice field 

4.37 

6 .44 

-0,117 

Vineyard 

4.48 

3.57 

-0.012 

Cereals 

7.57 

5.15 

-0,005 


For the calculations, radiometric value is coded with 8 bits (from oto 255) 
When studying the diagram, we can make following remarks ; 

- The projection of the rice-field class (1) on MSS5 axis will 
be merged with that of the vineyard class (2) : same density 
range in the chlorophyll band. 

- The projection on MSS7 3°ads to an ambiguity between class 1 
and class 3 (cereal) i existence of several irrigated cereal 
fields. 

This density slicing applied on MSS7 image will give rice field 
class mixed with cereal class. 

On the other hand, we can see that the plotted points which 
correspond to pairs of measurements x^, y^ in MSS5* MSS'7 bands lie along 
an axis. For each class, the slope of the axis depends on the correlation 
between x^, y^. If the characteristics were completely independent of one 
another, the mean value of their product would be 0, and the axes would 
be axes of x and y. 

In this example for class "2" and class "3", the MSS7 responses 
are rather independent of MSS5 responses. This means that the correlation 
coefficient, which is : 

2 x. y. 
i * i 

r = — * — is close to 0 

Cix.tfy 




64 
‘ * 
I”* 
84 I 

r 

i . 
i 
i 
i 
i 
i 
i 
i 

<54 I 


74 

<r 


84 

+ 


94 

• 4 * 


104 

* 


2 

2 


89 


2 

2 


2 2 


1 1 


33 1 
1 


2 2 2 " 

I 22 222 

I 222 2 2 
I 2 2222222 

I 222 22222 
99 I 2 2 22 22 

I 22 222 

I 2 2 22223 

I 222 2 

1 2 2 2 2 2 

104 I 2 2 22 

I 2 22222 
I 2222 2 
12 2 2 
I 222 22 2 1 

109 I 2223221 : 

I 22 3 

I 222 131 

I 2 2 21 

I 3 3 

114 X 33233 23 2 

13 3 3 

13 3 123 

I 


1 


1 


11 


1 


1 

31 

I 13 

II 111 
1311 11 

11 11 

21 11 11 
•si 1 T 

mi in 3ii 
31 3 11 ' 
131 1 1 

1 


1 ■ 
11 

1 

11 1 
11 1 
1 3 


1 


11 


33 


33 

3 

1 


3 1 


1 1 
11 3 li 
1 1 
1 11 


I 

119 I 
I 
.1 
I 
1 

124 1 
I 
I 
i 
i 
i 
i 
3 
I 
I 
I 
I 
I 
I 
I 

139 I 
I 
I 

' I* 


333 3 
333 
333 
3 
3 

333 


3 

333 

3 


3 33 

333 33 
3 


3 

3 


3 

3 


33 33 3 


129 


134 


33 33 
33 
3 
3 


3 
3 
3 
33 

33 3 
333 3 


3 33 




—i 

x 

-i 

X 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

X 

I 

I 

I 

I 

I 

I 

I 

I 

II 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

“I 


Figure 3 

Bispectral diagram of the responses of sampleS(MSS5 and MSS7 ban&s^ 
1 : ricefield, 2 : vineyard, 3 i cereals. 
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The slope of the axis can he calculated, if we assume that the cloud- formed 
fay plotted sample points of every class is an ellipse, by : 

0rt -2r.6x.(fy 

tg 26 “ — — ~ jL - 

ifx 2 - €y z 

with 0 : angle between x axis and the ellipse great axis. 

Since r is close to 0, either 0 will be very small, or else 
<5 x is nearly equal to<?y. 

In this case, where<Sx i (y y for the 2 classes, the possible 
conclusion is that 0 is very small, and the axes of class "2" and 
class H 3 ,r clouds nearly coincide with the x axis. 

On the other hand, when we compare the slopes of the axes of 
class "1" and class ,T 2”, as 

<5 xi ^<Sx 2 and 

<5 y\ >G y 2 : (dispersion of rice-field responses In MSS7 bandj 

since ri > r 2 ) 

tg28 i > tg2©2 

The slope of the rice field class axis is greatetf than that 
of the vineyard class. 

* 

Then, the fastest method of classification consist to choose 
a linear classifier in MSS5 _ MSS7 plane. The equation of the line can be 
easily determined on the diagram. 

The plane is then separated into 2 regions : region of the 
ricefield class and the region of other classes. We must however to ascertain 
that the region 1 does not include other classes (water bodies...). 


The calculations and notations are from = An Introduction to Mathematical 
Probability. Julian Lowell Coolidge - Dover publications, INC. NEW YORK. 
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.On other example of the application of the method is shown 
in the figure 4, where the plot of the bitemporal data permit to separate 
ricefield class from the others by a linear classification. 

XV - Result 

Figure 5 shows the result of classification using linear clas- 
sifier from the scene of July 6 1975. For 180 x 180 pixels (15 km x 15 km) 
this method gives 45 % area occupied by the ricefield. 

The calculation time is 0.07 sec CDC 7600. 

Figure 6 shows the results of classification method using 
Euclidian distance. On the same surface, the results give 41,3 % ricefields, 

20 % vineyards, 29 % areal and com. The cal .u la t ion time is 1,5 sec 
CDC 7600. 

For both results, the reference ricefields which are : Clos des 
Saules, Mas des Aubepines, Mas d’Adrien, Mas de Vert are well defined. 

In figure 7 the outline of a ricefield, the Clos des Saules, 
obtained by processing of LANDSAT imagery is compared with a ground survey 
map. The overall area of the field is found on the LAND5AT results, but in 
fact the roads and dikes within a block of parcels occupy a non productive 
area of from 15 to 20 % , according to investigations made on the 3 parcels : 
les Pebrieres, le Clos des Saules and le Mas de Vert. 

On the result of classification in figure 6, a smoothing 
technique is applied. For* every pixel to be classified, the 8 surrounding 
pixels are considered. The smoothing technique takes into account the distance 
between pixels, the class of the central pixel and that of the surrounding 
pixel s . 


This smoothing technique which gives good results in large homoge- 
nous zones is considered unsuitable for rice fields in the Camargue for 
identification of blocks of parcels. The figure 8 shows a degradation in the 
contour,* of the small area rice fields : this is indeed the case of the 
Mas des Saules, Mas des Aubepines, and the Mas d’ Adrien. 
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Figure 4 


Diagram of data from MSS7 band in 2 dates : Marcia 1973 and May 1973 
Class 1 : meadow and cereal, 2 : vineyard, 3 : Rhone, 4 : petit Rhone 
5 ; urban, 6 : ricefields 


*■••• 


mi ••• 








«« 

Mill 

UUII 

l«HM M 

Hnll 


IttttMM Mum <f 
•iflllill ••*•*••• 


till 

fill 


. 1C 

• » * 
14 1*1 
IMM 
IIU 


IM 


• I 


II 
lit 
11(11 
ttti r i 
i mi III 
Ml Mt 
I It 


«•« 

• •<1 
tllll 
Ml 


»•» 

Mil 


•tf *■««*•* 

tllliMM 

Ml 
• *C 
tilt 

• Ml 
« life 9 
I 


IMUMfl 
■ Ul 


• Ml III 

• f 

III III 

• < 1.1111 
• • Mt 
If 1*11 
life Mil 

at « 

it * • 


t« 

lit 


• M 

• Ml 

• Ml 

• MM* 

I ••II* 

limiM 


I MM 
III MM I 
• ••«< 
I • 


II II 

MIM 

• MM 

• MU 


Ml 

Ml 


Mill C 
ft «l 

I MMMI 
Ml Mil 


II 

• M I 

II II I 


II 

MM 

It 


1 1 1 4 1 

■ hum 

umiimm 


IM • 

M«tl I 


(Mil 

II *4*11* lltU*«Mt|MM 
I |MK if * %• I 


IMM 

«**• 

It* 


Ml 
■ •I* 

• * I • • 

• 4411 






imimmi 

IIMl'il 
ll\ M • Ml 
|IMM 
I'M 


l«U 

<muti 

min 

M*Mt 

• .till 

• lift 

Ml 

• ••I 


t I •«*••«»• 

llttlll'HflXI 

iimniMiim 

Ml I'M l*«t 

•»« MMMMMM* 


IM 

lilt 

IM 

IM 


i* • tinitMMiMt m*i<< • *• 


• •imiiui iMMtMiuu titi'tfM 

Mil 4M»MM I 114 MM( 

■ • «•!« II**# • **• 

t*M MitUml 

• ««« f IllMlfMI 


■ mu 

MM I 

M iix 

ItIMIf 

tM««M 


tfe 


IM 
• lit I 
M til 

If llll'lf Mil 

Mint i m”i tun 

VI* T* <411* M 1*111* l*l| 


it a 
imiim 
1 1 M * I 
1 ** 1 * 
t HIM 


• If 
•III 

• Mil 

• Ml 
MM 


V* 


M • «•■■• 

IM* II IM 
I* I'llill 
Ml IM MM I 
1 till 

IlMl Itttll 

IMM 1 1 * I • I 

fUlilllllllill 1 


**•1 


It 


IMMMlM 
IMWttll 
• II* MM 
III CM* 


It 


Mill 
• Mill 
Ml* 


• Mil 

•mi 

limt 

nmiit 

IIMIMMI 

iimiii 
| lMll.ll 
■ MMMtll 

• MMI 1* 
Ml Ml 

• •Mil* 

II *«**l 

III I l< 

III « 

Ml II 
IM 

l«l*»» 

1*1*11 
until 
1'1'Hl 
lunmi 
inm*« 
••*••••««’• mi 


• |M4* 

• 4 
HI 
MM 

• IMM 
t Mt Ml M 

1 •***«• 
Mil 


M*mMMI »I1 


M 
MM 
MM* 

• ••Ml 
Him* 
••••ft 
*•#♦«* 


Mill MIM 
l»*f KUI 
MM IMM 

im'iiii min is 

!•«««•• I 1 1 M 
•••It 
III I 
Mt t • 
Ml I 

• Mt 

• ••M I • 

• MMMI f 


• •••Ml I I 

IMIMMI MM 

Mill 


Ml 

«•« •■••«• 

■ Ml • Ml 




4 • • * • •• It 


M III! 
I* till 


• Mf t 
• lllil 
I • «#*• 

If It* Ml*** 

• Ml II • 

■ Ml •• 

Mil t 

It MMI 

| it 1 1 life 

| 1 1 1 • ■ 
tv* tmi 

llifrmmim » 

t *«••**« mm i« mint 

• ««««•»«* *«•««•• m.»i#»n i»»#*t#*** , »" 

•imi«*i iM«t •i»m m ,.ii *“''* *!!!!!. 

■*••1* ««•««••« •«•«•* i* **•••••« •* •*•••*•#** 


•««*• 
M« Ml 
mi i 

lt«l 

IM 

111 


*(•• 
M tf • 
«•• 






* *«••• 

• •««•! 

* t(ll« 

•iimiiiu 


IM 

• Ml 

• It 
Ml 


• I • 


•••It* 

mmmi 

• t c •••tiiii 


Figure 5 


Classification of July 1975 scene using linear classifier 
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Classification of the July 6 1975 scene using Euclidean distance 
dark : rice 

medium : vineyard 
light : cereal 
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Figure 8 

Result from classification method on July 6 1976 imagery followed 

hy smoothing technique 
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Unsupervised classification method applied to rice-field 
characterisation from LAND SAT 2 imagery - " 

P, CASSIBAME, T. LE TOM 
Centre d ’Etude Spatialo des llayonnements 


Abstract 

The mobile centers clustering technique used by the C.1S.S.R. as unsupervised 
classification, method is described and especially the initialisation of 
the process which is carried out in such away to reduce the processing time 
and the memory needed, is explained. 

The application of the method for ricefield mapping from LANDSAT data is 
shown, and the results are compared with that obtained by supervised techniques 
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The evolving of supervised methods is based on the existence and 
the properties of prototype elements -which correspond to ground objects, 
determined by ground investigations. 


"Whenever prototype samples are not acquired, unsupervised methods 
must be developed. 


The technique evolved by the C.E.S.R. is a "mobile centers" cluster- 
ing technique called "Nuees Dyn antiques" (Dynamic Clusters) by this author, 
E. DIDAY (1) . 


I - DESCRIPTION OF THE METHOD 


I . 1 . Principle 


Given a metric space (E,d) and a set k £ J called the set 
of centres. 


To k is associated a partition of E into polyedral, convexe subsets 


Ec : 


£c = 1 ze £ / Ve'^e £ k : d^,e) ^ d 


Each set Ec has a center of gravity Gc and for each Ec, a cluster made 
up by points Nc called nuclei is built. Nc set is composed by points condens- 
ed aroung Ge and is determined by its cardinal nc : 


We = \ St' £ E, 


Wc = <hc , V fc 4- Nc ' d (ct,G t 

4 

The center of gravity of Nc cluster is then calculated ^ Gc and another 
centre K .1 is defined. The e 
around the new centre of gravity G 


i !■: 


i r 


if. 


\ i 

: f 


<4 


set of centre Ki is defined. The elements of which are the nearest elements 

, 1 * 

V « 

C '* 


A new partition of E into subset 1 E 1 -center G*~ is then constructed 

c 1 c 

as above end the subset Gc is then superseded by G c * 
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The dynamics cluster method consists then to iterate the process, 
i.e. to build new centers Gc, G^, G c ..., G c and to stop at the iteration i 
if g£ - G * +1 \/c. 

The method is justifies by the fact that the passage from G c to G^, 

2 

G c , ... reduces the intra-classes variance of Ec. 

1.2. Convergence of the process 

If the cluster E is finite, there is then a finite number of partitions 
of E. To each partition corresponds a total internal variance that decreases 
strictly after every iteration. 

The variance will reach a minimum after a certain number of iterations. 
Further, successive iterations modify only boundary points, i.e. the when the 
distances from these points to 2 distinct centers are equal. 

Remark 1 : Instead of the center of gravity G c of N c set an element 
of N c set can be taken. The creation of any artificial non convexe class is 
then avoided. 

Remark 2 : E. DIDAY advises to use all the points of the N c set to 
build the new partition E* . That allows : 

. to reduce the effect of the marginal elements. 

. to apply various metrix, for example the Mahalanobis metrix which 
cannot be used when the N c set is reduced to one point. 

. to avoid the clustering of the 2 points far from each other but 
linked by a line of points. 

. to use the totality of the samples if the method is used in supervised 
manner. 

1.3 . Initialisation of the process 

The determination of the set of centers, k, can be made : 

- either in studying the results from a correspondence analysis and 
that needs a supplementary treatment. 

- or randomly, and that leads to a great number of iterations 

RF- - 'J‘- 1 - ^ 

— > i. 
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To avoid these 2 inconvenient s, the technique carried out by the 
C.E.S.R. is the following : 

The area under study is subdivised into sub-areas Z lt Z 2 Z n , 

each of which is a square of 60 x 60 pixels. 

The process is applied first to the zone Z\. The initial set of 
centers is chosen randomly and the iteration number is limited to 5. 

If after 5 iterations, the convergence of the process is not obtained, 
it is decided that there was bad initialisation, then an other random 
choice of initial centers is made and the process starts again. 

When the convergence on Zi is obtained, sone Zi+1 is then studied 
with the last calculated set of centers of Zi as initial centers. As 
there is close neighbourhood between Zi and Zi+1 with best chance of 
fitting and the convergence is reached practically after one iteration. 


II - APPLICATION TO AGRICULTURAL MAPPING IN CAMARGUE 


The test-site is an area of 180 x 180 pixels (about 15 km x 15 km) 
in the region surrounding Arles (figure 1), where the ground truths are 
well known by ground investigations. The results from unsupervised class- 
ification can be then compared with that obtained by supervised method. 

The example shown in the following concerns multitemporal Landsat 
data. The MSS5 and MSS7, bands from the two scenes of March 21 1S73 and 
May 14 1973 are used. 

II . 1 . Ground truths 

The main classes existing in the area are : 

- The Rhone 

- The Petit Rhone 

- ricefields 

- vineyards 

meadows and cereals 

- Arles 


• * • / . . . 





At the date of March 21, the rice-fields were ploughed, wet parcels, 
the vineyards were at the beginning of the vegetative state and the fields 
of cereals as well as meadows, formed a green carpet. However, other fields 
could look like ricefields (ploughed com fields). 

At may 14, ricefields are flooded, except some fields where the flood- 
ing operation was done on May 15. 

II. 2. Application of the unsupervised method 

The area under study consists of 60 x 60 pixels subareas. 


The process is starded on Z^, where 6 initial centers 
are picked randomly. To the 3600 points of Z ^ are applied 
the clustering algorithm with a chosen metrix. 


The metrix used in this example is the euclidean distance. 


*11 

Z 12 

Z 13 

Z 21 

Z 22 

Z 23 

Z 3! 

Z 32 

Z 33 


When the convergence is obtained in Z^, the results of the classifica- 
tion for 60 x 60 points of Z ^ is preserved for display, and the last 
computed centers are then conserved to be taken ad initial centers for Z ^ 

In the case of LANDSAT imagery, each sub-area corresponds to a zone of 
about 4,8 km x 4,8 km. Thus, geographically, two contiguous subareas are in 
most case included in a same "ecologic region". This means that the same 
types of vegetation, of land-use exist very probably in two contiguous sub- 
areas . 


Practically, one or two iteration is necessary for computing Z ^ with 
the last centers of Z ^ as initial centers. It can occur that "landscape" 
changes suddenly -when the sub-area is on the boundary of 2 "ecologic regions" 
Then other itererations are needed, or an choice or random initial centers 
must be done. 
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The contiguity notion can be applied in both directions line and 
column (geographically North South or East-West). Thus, for the 

initial centers came from Z^.,. 

Figure (2) shows the results of unsupervised classification obtained 
by display on Comtal device of the 9 subareas. The six classes, represented 

by 6 colours, are interpreted as follows : 

- Class 1 : dark blue : rice field 

2 : pink : meadow and low state cereals 

3 : light blue : ? 

4 : yellow : Rhone 

5 : green : vineyards 

6 ; red : petit Rhone and/or Arles ? 



This map is compared then with that obtained by a supervised method 
using the same distance (Euclidean). 

On figure (3), the following classes are taken into account, after 
a study of reference samples : 

- Class 1 : dark blue : Rhone 

2 : pink : petit Rhone 

3 : light blue : ricefields 

4 : yellow : vineyards 

5 : green : meadow and low cereals 

6 : red : Arles 

7 black : unclassified points 

The computed map obtained by supervised method is studied and compared 
with ground truths. ^ 

Classes 1, 3, 4, 5 are in good concordance with the ground truths, 
especially rice-fields which are delimited. Class 2 (petit Rhone) is partly 
mixed with class (Rhone) and class 6 does not represent the urban area (Arles) 
as expected. 


. . . / . . . 
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The comparison of the 2 computed maps ledd then to the following 
remarks concerning unsupervised computed map (fig 3). 

- the ricefields are well placed but their contours are not well 
delimited, especially in the North West side of the test-site. 

The edges of the Rhone are partly computed as rice-fields. 

- Meadows and cereals have like ley the same delimitation in the 2 maps. 

- The vineyard class does not correspond to the totality of the vine- 
yards in the areas and class 3 is not identified, it should corres- 
pond to an intermediate radiometric responses cluster between vine- 
yards and meadows and cereals. 

- The urban area, in both maps, cannot be isolated. 

II . 4 . Processing time problem 

The processing time, on the base of CP time CDC 7600, is 0,9 sec for 
the 180 x 180 pixels test-site with supervised method using Euclidean distance. 

For unsupervised method, the processing time, which depends on : the 
applied metrix, on the number of elements of every nucleus -and on the number of 
iteration, is evaluated in this particular case as follows : 

For nucleii of 10 elements, the processing time for one iteration is 
0,25 sec CDC 7600 for the first subarea -the metrix used is the Euclidean one. 

The processing time for 9 x 60 x 60 pixels is approximately : 

0,25 x 5 + 0,25 x 2 x 8 = 5,25 sec 



5 iterations 2 iterations (average number) : 

- when Mahalanobis or X 2 metrix is used, the processing time is to 
multiply by 5/4. 
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Figure 2 

Computed map obtained by "mobile centers clustering technique" 
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Figure 3 

Computed map obtained by supervised method 



